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Signals of dark matter in a supersymmetric two dark 

matter model 



Abstract 

Supersymmetric radiative neutrino mass models have often two dark matter candidates. 
One is the usual lightest neutralino with odd R parity and the other is a new neutral par- 
ticle whose stability is guaranteed by a discrete symmetry that forbids tree-level neutrino 
Yukawa couplings. If their relic abundance is comparable, dark matter phenomenology 
can be largely different from the minimal supersymmetric standard model (MSSM). We 
study this in a supersymmetric radiative neutrino mass model with the conserved R par- 
ity and a Z2 symmetry weakly broken by the anomaly effect. The second dark matter 
with odd parity of this new Z2 is metastable and decays to the neutralino dark matter. 
Charged particles and photons associated to this decay can cause the deviation from the 
expected background of the cosmic rays. Direct search of the neutralino dark matter is 
also expected to show different features from the MSSM since the relic abundance is not 
composed of the neutralino dark matter only. We discuss the nature of dark matter in 
this model by analyzing these signals quantitatively. 
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1 Introduction 



The explanation of small neutrino masses [T] and dark matter [2j seems to be a key 
ingredient to consider physics beyond the standard model (SM). An interesting possibility 
of such extensions may be models which can closely relate neutrino masses to dark matter 
(DM). In this kind of models, a discrete symmetry is often introduced to forbid tree level 
Dirac neutrino masses. Some of additional particles introduced to commit neutrino mass 
generation have its charge such that it can forbid the lightest one to decay into the SM 
particles. This stable particle becomes DM. This DM is a crucial ingredient of the neutrino 
mass generation in this scenario. 

The radiative seesaw model proposed by Ma [3j is its simple and interesting exam- 
ple]^ Both the numbers of new particles and free parameters are comparably small. Its 
supersymmetric extension is also straightforward [IHl [H] Ifl Moreover, if we introduce an 
anomalous U(l) symmetry in this extension [13], we could explain the origin of the dis- 
crete symmetry, required hierarchical structure of both couplings and masses due to the 
Frogatt-Nielsen mechanism [T6] . Both neutrino oscillation data and DM relic abun- 
dance can also be explained consistently with lepton flavor violating processes such as 
fi — )■ 67. A characteristic feature in such an extension with R parity conservation is that 
the model has two DM candidateqj. One is the lightest superparticle whose stability is 
guaranteed by the R parity. The other one is a new particle introduced for the neutrino 
mass generation and its stability is guaranteed by the new Z2 symmetry. As a result, the 
model shows discriminative differences from the ordinary minimal supersymmetric SM 
(MSSM) in the DM search. For example, if the recently reported cosmic ray anomalies 
pn [T8] are considered as the DM signature of the model, they may be explained not by 
the DM annihilation |19l |20l |2T] as in the MSSM but by the DM decay [221 [23l |Ml |25l [26] . 
In fact, if the Z2 symmetry is violated by the anomaly effect, the DM guaranteed its sta- 
bility by the Z2 symmetry can decay into the lightest neutralino [TH [I3| . Direct search 
of the DM could also show the difference from the MSSM. 

In this paper, we study signals of the DM in the supersymmetric extension of the 



lot of radiative neutrino mass models exist now. Phenomenology including the DM nature in such 

models has been studied in a lot of works [H O [SI [71 [H [3] ■ 

relevant supersymmetric model is also considered in a different context in |12j . 
^Multicomponent DM and its phenomenology are studied in a different model [15] . 
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Table 1 Matter contents and their quantum number. Z2 is a remnant symmetry of the assumed 
anomalous U(l) caused by the symmetry breaking at a high energy region. 

Ma model. The model is considered as an effective model due to spontaneously broken 
anomalous U(l) gauge symmetry. It naturally brings the weakly broken Z2 symmetry to 
the model in addition to the conserved R parity. We discuss signatures due to the decay 
of the unstable DM and also the direct detection of the DM through the elastic scattering 
with nuclei. 

The paper is organized as follows. In section 2 we address the model and explain the 
nature of the DM sector which is imposed by various experimental results. In section 
3 several signals expected in the DM sector are analyzed. In particular, the decaying 
DM is studied to explain the cosmic ray anomalies reported recently. A feature of the 
monochromatic gamma yielded through the DM radiative decay is also studied. Finally, 
we discuss the direct search of the DM. Section 4 is devoted to the summary. 

2 A supersymmetric radiative neutrino mass model 

The radiative seesaw model proposed in [3] is an extension of the SM with three right- 
handed neutrinos and an inert doublet scalar. The latter is assumed to have no vacuum 
expectation value (VEV) and no coupling with quarks. Although the model is very simple 
and has several interesting features [H IS |6] , it has some faults, that is, the existence of an 
extremely small coupling and the ordinary hierarchy problem. These may be improved 
by extending the model with supersymmetry and an anomalous U(l) symmetry \13l. We 
focus our present study on this model, which has a Z2 symmetry as a remnant subgroup 
after the spontaneous symmetry breaking of this anomalous U(l). Matter contents of the 
model and their Z2 charge are summarized in Table 1. 
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The most general superpotential invariant under the imposed symmetry is 
W = hYjQiU:;Hu + h^^Q,D''^Hd + hfL,E^Ha + fiHHuHd 



(1) 



+ QMpie ^'LiTju. 

This can be obtained as the low energy effective theory through the spontaneous breaking 
of the anomalous U(l) as shown in ^13j. The last term in W is induced by an anomaly 
effect [SniET]. This term breaks the Z2 symmetry very weakly if bi is large enough. Since 
the Z2 symmetry is not exactly conserved, the lightest field with odd parity of the Z2 is 
unstable. However, the lifetime can be longer than the age of universe and it behaves 
as the DM. Thus, we have two DM components in the model as long as the R parity is 
conserved. 

Soft supersymmetry breaking terms associated with the superpotential W are intro- 
duced as follows, 



C 



SB 



+A(/ipiiV,% + XuVuHd4> + XdfjdHu4> + h.c.) 



B UnVuVd + 



+ -MiN^' + CiMpie-'^Lif]^ + h.c. 



(2) 



The scalar components are represented by putting a tilde on the character of the cor- 
responding chiral superfield except for the ordinary Higgs chiral superfields and if^. 
Universality of soft supersymmetry breaking A and B parameters is assumed, for simplic- 
ity. Moreover, we confine our following consideration to the case where soft masses for 
all the scalar partners are flavor diagonal and universal unless we mention it. They are 
denoted by uiq. 

Neutrino masses are generated through the one-loop diagram as discussed in [13]. If 
we focus our attention to the special flavor structure for neutrino Yukawa couplings such 
as [5] 
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the neutrino mass matrix is found to be expressed as 
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(4) 



= 00^2' [di^i^ - 9{Mu rrin. 



This mass matrix induces the tri-bimaximal MNS matrix^ Mass scales for the neutrino 
masses are determined by Aj, which is defined as 

ml - ml + ml ln(mg/mg) - _ A»Arftan/j 
{ml-miy ' ~ l + tan2/3' 

where = f sin /3 and {H^) = v cos (3. Xu,d are assumed to be real, for simplicity. m^_|_ 
are the mass eigenvalues of the neutral scalar components of 1]^^^, which are defined as 
m^_i_ ~ /i^ + ml ± Bfijj. If Mi and m^-t have the values of 0(1) TeV, mass eigenvalues 
of neutrinos can be suitable values as long as Xu and A^ take very small values such as 
A„A, = O(10-«). 

Before proceeding the analysis of the DM phenomenology, it is useful to address free 
parameters in the neutrino sector of the model. The relevant parameters are summarized 
as A, firi, mo, B and also \hi\, ipi, Mi {i = 1,2,3). We restrict our study to the case 
with Ml ^ M2 < M3 which allows the coannihilation of ipN-i and 11)^2 (the fermionic 
components of A'^^^ and N2). We consider this case since it brings an interesting aspect 
in DM phenomenology as seen below. Since one eigenvalue of (jlj) is zero, the neutrino 
oscillation data tell us that remaining eigenvalues should be \/ Am^^j^ and y^Am^. This 
imposes the parameters to satisfy the relations 

Thus, after using these relations, the free parameters in the neutrino sector can be confined 
to 

A, Ml, M3, /i^, mo, B, ifi. (7) 

Here we search parameter regions consistent with the experimental data for the lepton 
sector. For this purpose, we can use neutrino oscillation data [I] and the constraints 
from lepton fiavor violating processes (LFV) such as Br(/i — 67) < 1.2 x 10~^^ [27] and 
Br(r —7- yU7) < 4.4 x 10~* [25]. We fix a part of the parameters listed in eq. ([7]) aq^ 

A = 1.16 X 10-^ mo = 400 GeV, M3 = 9000 GeV, - = 0. (8) 



^The charged lepton mass matrix is assumed to be diagonal when we consider the flavor structure of 

neutrino Yukawa couplings ([3]). 

* In this analysis we use the parameters different from the ones used in [T3]. It could cause some 

differences for the bounds of parameters between these two cases. For example, since we use a larger value 
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Fig. 1 The left frame shows a parameter region in the (Afi,/i^) plane which is consistent with the 
neutrino oscillation data and the LFV constraints for B = 250 GeV. An upper red solid line represent a 
contour for Br(/i — > 67) = 1.2 x 10^^^ which is the present experimental upper bound [27 and another red 
solid line represents a reference value Br(/i ej) = 0.6 x 10^^^. Blue dotted lines represent the contours 
of Br(T fi-f) for reference values 2 x 10^^ (the left one) and 0.6 x 10^^ (the right one). Its experimental 
upper bound is 4.4 x 10^^ 28J. Mi = is satisfied on a black dashed line. In the right frame these 
branching ratios are plotted as a function of B for AIi — 3200 GeV and /i,, ~ 3600 GeV. It shows an 
allowed region of B which is consistent with the neutrino oscillation data and the LFV constraints. A red 
solid line and a blue dotted line represent Br(/i — >■ 67) and Br(r — > /i7) in this model, respectively. Thin 
black dotted lines represent the experimental bounds for Br(// — > ej) (the lower one) and Br(T — > jij) 
(the upper one). 

If we use these parameters in eq. (|6]) and the formulas for the branching ratio of the LFV 
jl3] . we can find a parameter region in the (Mi,/i^) plane which is consistent with both 
the neutrino oscillation data and the constraints from the LFV. It is plotted in the left 
frame of Fig. 1 for B = 250 GeV. In this figure, a region sandwiched by the upper red 
solid line which represents the contour Br(/i — )■ 67) = 1.2 x 10~^^ and the black dashed 
line which represents Mi = is an allowed region if tp^i is assumed to be the lightest 
Z2 odd particle. It shows that the LFV constraints can be satisfied for /i^ ^ 4850 GeV. 
This result does not depend on the mo value sensitively in the region where Mi, /i^ ^ mo 
is satisfied as long as B is fixed in the region B ^ Mi, yU^. If we take a larger A, the 

of A/a than the one in p3] here, the LFV constraints are satisfied for a smaller mg value compared with 
the one discussed there. These parameters are adopted here since they are favorable for the explanation 
of the cosmic ray anomalies as discussed later. 
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condition (E]) can be satisfied by smaller neutrino Yukawa couplings. In that case we note 
that the LFV constraints become weaker. It is interesting that this allowed region in the 
{Ml, fin) plane relevant to the following analysis is within the reach of /i — )• 67 search in 
the MEG experiment. It aims to search for /x"*" — )■ e"*"7 decay with a sensitivity of a few 

In the right frame of Fig. 1, we plot the branching ratio for yU — > 67 and r — )■ /i7 as a 
function of B for Mi = 3200 GeV and fi^ = 3600 GeV which are contained in the allowed 
region as shown in the left frame. This figure shows that the LFV bounds can be satisfied 
for B ^ 200 GeV. If we fix B to 250 GeV for example, Yukawa couplings are found to 
have rather large values such as \h1 + h"^]^^"^ ~ 2.98 and j/isl — 1.14. Although the values 
of Yukawa couplings gradually decreases for larger values of B, they are always large in 
this figure. These large Yukawa couplings are required to reduce the relic abundance 
of V'TVi with such a large mass sufficiently]^ Although they could cause a problem for 
perturbativity of the model, we can escape this fault of the model by considering the 
phases of neutrino Yukawa couplings. This point is discussed below. 



3 Signals of the DM 

3.1 Relic abundance of two DM 

The model has two types of DM candidate. One of them is the lightest neutralino x 
whose stability is guaranteed by the R parity as in the case of the MSSM. The other one 
is the lightest neutral state composed of the components of Z2 odd chiral supermultiplets 
N^, rj^^ and (p. In the following study, we assume that tpNi (the fermionic component of 
N^) is the lightest one among these candidates. Since this Z2 is not an exact symmetry 
but is weakly broken by the last term of W through anomaly effect, ipNi is not stable. 
However, it could have a longer lifetime than the age of the universe as long as bi is large 
enough. If this is the case, the DM relic abundance suggested by the WMAP [2] should 
be satisfied by both of these contributions. This condition is expressed as 

n^h"^ + n^^y = 0.11. (9) 

^If we assume smaller values for Mi and /x^, small neutrino Yukawa couplings can explain the ipNi 
relic abundance consistently with other constraints by fixing A to a larger value. 
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Fig. 2 The relic abundance /i^/O.ll as a function of the mass of i/jn-^ for the cases B ~ 

200,250,300 GeV. (fi — ip2 is fixed to (the left frame) and j (the right frame), respectively. A are 
fixed to 1.16 X 10"^ and 7.74 x 10^^ in the left and right frame, respectively. Other parameters used to 
draw these figures are explained in the text. 



If fi^ ^ ^ipN-^ is satisfied, DM physics is the same as the one of the MSSM. However, we 
are interested in a different situation from the MSSM, where both of them cause the same 
order contributions. In order to study DM physics for such a case, we search a parameter 
region which brings this situation within the parameter space discussed in the previous 
part. 

First, we consider the annihilation processes which determine the relic abundance of 
ipNi- The annihilation is induced through the t- and u-channel r/u exchange. If tpN2 has 
the almost degenerate mass with tpNi, we need take account of the coannihilation effect 
[32] . We suppose such a situation here. The possible final states of such processes are 
composed of a pair of lepton and antilepton or a pair of slepton and antislepton. Applying 
the method developed in [32l [33] to this model, we can estimate the relic abundance 
Q^^^h"^. The details can be found in [13]. In this estimation, we use the parameters given 
in eq. ([8]) and /i^ = 3600 GeV which can be consistent with the neutrino oscillation data 



and the LFV constraints as seen before 
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In the left frame of Fig. 2, the relic abundance f^^'jvi plotted as a function of the 

ipN-^^ mass Ml for typical values of 5. In the right frame we also plot the same figure for 
A = 7.74 X 10^^ in the case of (pi — ip2 = j- Other parameters are fixed to the same values 



^'•^ Since the difference between Mi and is 10% in case of Mi = 3200 GeV, jtiq — 400 GeV and 
B = 200 GeV for example, the coannihilation of -0^1 and ry_ could play some role for a larger B [35]. 
However, we neglect their coannihilation effect in this analysis. 
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Fig. 3 Contours of the mass of the hghtest neutrahno x (left frame) and contours of its composition Zu 
(right frame) in the {fiH,Mg) plane. Only a region with fin > 500 GeV is allowed for mo = 400 GeV 
and B > 200 GeV. 

as the ones in the left frame. As noted above, since A is fixed to the larger value compared 
with the one used in the left frame, neutrino Yukawa couplings can take smaller values 
keeping the consistency with the condition ([6]). On the other hand, ipi ^ generates the 
s-wave contribution in the coannihilation cross section for ipj^^^ and t/'^Vj [13]. As a result, 
neutrino Yukawa couplings such as \h\ + hW^^"^ ~ 1.19 and |/i3| ^ 0.45 are sufficient to 
realize the suitable relic abundance in case of Mi = 3200 GeV, for example. These values 
of Yukawa couplings are much smaller than the ones in the ipi = <^2 case. Because of 
this feature, the present LFV limits give no constraints on the model in this case. This 
situation is largely different from the ipi = ip2 case where the LFV constraints could play 
a crucial role to restrict the parameter space as shown in Fig. 1. 

Next, we consider the relic abundance of the lightest neutralino which is defined by 



X 



(10) 



where Zu (z = 1 ~ 4) are determined by diagonalizing the neutralino mass matrix 

/ — cos /3 sin Owniz sin /3 sin Ow^nz ^ 

- sin (3 cos 9]v^z 



M, 



N 



My^ COS /3 cos Ow^z 

— cos (3 sin 6wmz cos (3 cos Owmz — /i// 

y sin/3 sin — sin /3 cos 6'i^m2 — yU// j 

The annihilation of x occurs through various processes depending on its composition Zu 
in the same way as the MSSM. Final states of the % annihilation are composed of all 
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Fig. 4 Contours of the relic abundance H^h^ 



in the {p,H,Mg) plane. Only a region with 



Hh > 500 GeV is allowed for mo = 400 GeV and B > 200 GeV. 



the SM particles which are hghter than x- The favorable parameter regions to explain 
the DM abundance by x have been studied in detail [31]. We follow such studies in the 
present case keeping in mind the consistency with the parameters used in the neutrino 
sector. 

The relic abundance of x is determined by the mass and the composition Zu which fixes 
the interaction of x with the SM particles. The parameters relevant to them are tan /3, fin, 
soft supersymmetry breaking parameters rriQ, A, B and the gaugino masses M^y, M^. Here 
we note that some of these parameters are related each other. M^^ = 2M^ is expected at 
the weak scale from the unification relation among gaugino masses. Since B is required to 
satisfy an electroweak symmetry breaking condition B = {th?^^ +mjj^ + 2/i|^) sin 2/3/2/i//, 
it is determined by fin if we fix the values of rriH^, tan/3. In this analysis rriH^ 

and rriHa are fixed to 500 GeV. We also take tan /3 = 10 which can be consistent with 
the discussion on the neutrino sector by tuning the value of A^A^. Since other parameters 
are fixed to proper values at the weak scale, fin and are treated as free parameters. 
Numerical calculation is executed by using the public code micrOMEGAs [35J. 

In Fig. 3, we plot contours of the mass and the composition Zu of x i^i the 
{fiH, M^) plane. When we see this figure, we have to remind that the allowed region 
should satisfy a condition < mo, which is required since x is DM. Here we consider 
a case with mo = 400 GeV and B > 200 GeV which are used in the study of neutrino 
sector for Mi = 3200 GeV. In this case fin > 500 GeV is required by the electroweak 
symmetry breaking condition given above. Thus, Fig. 3 shows that x is bino dominated 
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at the region with < 400 GeV where x can be the hghtest superparticlei I In Fig. 4, 
the contours of the rehc abundance Q^h"^ of x is plotted in the {fin, Mj^) plane. As found 
from Figs. 3 and 4, x can be a DM component with the substantial abundance at the 
above mentioned bino dominated region. Since the ip^-^ abundance is not sensitive to the 
value of mo if rriQ < is satisfied, we can find parameters for which x is a DM component 
with substantial abundance under the condition ([9]). It is useful to note that in this model 
X could be an important DM component in the region where it is rejected as the DM in 
the MSSM framework. 

3.2 Decay of the right-handed neutrino dark matter 

In the previous part we showed that the DM can be composed of two components ■^/'at-^ 
and X which have the same order abundance. However, ip^^ is not stable since the Z2 
symmetry which guarantees its stability is not exact. Since this symmetry is considered 
to be a remnant symmetry left after the spontaneous breaking of the anomalous U(l) 
at a high energy region, it is broken by the anomaly effect. In fact, the Green-Schwarz 
anomaly cancellation mechanism induces the Z2 violating interaction as the last term of 
W nonperturbatively. If ■j/'tvi is heavier than x, this interaction brings the decay of ■j/'at^ 
to X through the diagrams shown in Fig. 5. 

In order to examine whether ipi^^ can be dealt as the DM, we estimate the lifetime of 
'j/'TVi due to the decay caused by this interaction. It can be roughly estimated as 

where we use \hi\, Cj ~ 1 and Mi ^ ttiq. From this formula, we find that ipi^^ can have 
a sufficiently long lifetime compared with the age of the universe, as long as hi > 79 is 
satisfied. Thus, although the true stable DM is the lightest neutralino we need to take 
account of the contribution of tpNi to the relic DM abundance and investigate the DM 
phenomenology. 

Charged particle observation in the cosmic rays by PAMELA [T7] and Fermi-LAT [T8] 

suggests that there are deviations from the expected background. The possibility has been 

discussed that these are consequences of the DM physics. However, if we consider that 

^^In the Mg > 400 GeV region, DM is a sneutrino which is difficult to realize the right relic abundance 
because of its effeicient annihilation. 
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Fig. 5 Decay processes of i/jatj to x- A bulb represents the anomaly induced interaction CiB Mp\e~''' Lifju- 

they are yielded by the annihilation of the DM, the annihilation cross section required 
for the explanation of the relic abundance is too small [19j. Some enhancement of the 
annihilation cross section at the present universe seems to be necessary [201 EI]- On the 
other hand, if we consider the decay of the DM, these anomalies found in the cosmic rays 
can be understood as long as its lifetime is sufficiently long [221 123] . 

In the present model, particles yielded in the decay of ip^-^ may bring the required 
extra contributions to the cosmic rays. The expected flux depends on the scale of the 
Z2 breaking CiBMp\e~^^ in eq. 02]) [111 [IS]- Iii fact, if 6, ~ 88 is satisfied, the anomaly 
induced interaction causes a large enhancement factor of 0(10^^) in eq. f lT2|) to realize 
a long lifetime of 0(10^^) sec for ipNi- This lifetime is known to be suitable to explain 
the anomalies found in the charged cosmic rays discussed above. Moreover, since ipisf-^ 
couples only with leptons and sleptons, this decay could yield only leptons and photon 
other than x- The fiavor structure of neutrino Yukawa couplings ([3D can restrict the final 
charged leptons to ^ and r dominant lyi^i This feature makes the model favorable for 
the explanation of the above mentioned cosmic ray anomalies. In the following parts, we 
study the positron and electron fiux predicted by this decay process and compare it with 
the data. We also estimate the nature of photon fiux expected in the radiative decay of 
ipNi which is shown in Fig. 5. 

The metastable ipNi decays to x^a^/3 through the left-handed diagram in Fig. 5. This 
decay is caused by the anomaly induced ia-V mixing. Since the mass of ipNi is of 0(1) TeV, 
M|(^ mo) < Ml is naturally expected from a view point of low energy supersymmetry. In 
this case the intermediate slepton £a is considered to be produced as an on-shell state. To 
take account of this possibility, we use the propagator for the slepton which contains 
the effect of decay width . 



^^Since x is dominated by the bino in our considering parameter region, we need to impose Cg ^ in 
the anomaly induced interaction CiBMp\e^^' Lirju additionally in order to keep this feature. 
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The differential decay width to the final state ia is given by 



dV 



deF{e)- 



2EaEi3 
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(13) 



dE^'"^"' . 4(47r)3Mi Jo '^"^ '"'Ml + £;«(cose- 1) 

where the 6 is the angle between ia and ijs. This formula is symmetric for the exchange 



of ia and i/s. The spin averaged amplitude |A^a/3P is expressed as 
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In these formulas we use the definitions such as 



(14) 



F{e) = sin^ + ( 



-K — 9) cos9, 



E, 



Ml -ml- 2MiEa 
2[Mi + Ea{cos9 -1)] 



E^ 



El + El + 2EaEBCOs9 + ml 



X' 



t = Ml- 2MiEa, u 



Ml - 2MiEp 



(15) 



where g' and g are the gauge coupling constants for U(l)y and SU(2)j^, respectively. The 
mixing angle 6'^ between r]u and r^^, can be taken as = 7r/4 since their soft scalar masses 
are assumed to be universal In this derivation we use the universality for the slepton 
masses and 6,. We also assume the fiavor independent slepton decay width F^. Since we 
use the fiavor structure ([3]) for neutrino Yukawa couplings and Cg = 0, sufiices a and /3 
in eq. (fT3il run over the lepton fiavor and r. Thus, the decay of ipN^ does not yield 
positron directly in the final state. 

In this model the positron is generated through the decay of yU"*" and r^. In the 
following positron fiux calculation, we use the positron spectrum '^^^'j^^ obtained from 
the simulation by using the MONTE CARLO code in the public package PYTHIA [36] which 



Even if they are not universal, however, this is a good approximation as long as /i^ is much larger 
than their soft scalar masses. 
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can generate the lepton whose energy distribution is given by eq. ( 1T3|) and calculate 
the positron from the decay of this lepton. The positron spectrum obtained through this 
calculation is shown in Fig. 6. Although two leptons are contained in the final state of the 
ipNi decay , the final positron flux is estimated by summing the contribution from each 
lepton which can be treated independently based on eq. (IT^ . Thus, using this spectrum, 
the positron flux yielded through the decay of ipNi is expected to be observed at the earth 
as 122] 



dE'Ge+iE,E') 



Br. 



dE' 



(16) 



where E^ax = (Mf - m2)/2Mi and Br^^ = Y^e^ ^(''pNi L^(3X)/'^tot where Ttot is the 
total decay width of ipNi, which is fixed by including the final states with neutrinos. This 
Br^^ is almost determined by h^^ and Cq, since we assume that the slepton masses and ba 
in the last term of the superpotential W are universal. Moreover, Br^^ is determined by 
Ca only in the present case since the flavor structure (E]) for h^^ is adopted. Values of 
and Br^^ used here are shown in Table 2. 

The positron Green's function Ge+ can be approximately written as 



V'iVi 



10 



16 



exp[a + b{E^-^ - E"-^)] cm^^ s 



E^ 



(17) 



where a, b and 6 depend on the diffusion model and the assumed halo profile [221 EHl [39] . 
Since the result is known not to be heavily dependent on these in the case of DM decay, 
we use the MED model [39j and the NFW profile [lO]. They fix these parameters to 
a = —1.0203,6 = —1.4493 and 6 = 0.70. We also assume that the two DM components 
ipNi and X have the same density profile in our galaxy. 





case (a) 


case (b) 


case (c) 


Ce • . Cj- 


0:1:1 


0:1:0 


0:0:1 


Bre+ : Br^+ : Btt-+ : Br^— 


■ 1 -i • i 

^ • 4 '4 • 4 


" • 12 • 12 • 12 


^ • 12 • 12 • 12 



Table 2 Values of Cq and Br^^ used in the calculation for the positron. It should be noted that Br^4 
(Br^±) takes a nonzero value even if (c^-) is zero. 
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Fig. 6 The positron spectrum 



calculated from ^ by using PYTHIA. Mi = 3.2 TeV and 



(IE 



300 GeV are taken as a typical example here. 



The background flux for electrons and positrons in the cosmic rays are given by 

0.16E-1-1 



prim.bkg 



sec.bkc 



1 + 11^0-9 + 3.2^2-15' 

0.7E°-^ 

1 + 110^1-5 + 600^2-9 + 580^4-2' 
4.5^0-7 



1 + 650^2-3 + 1500^4-2 ^ ' 

in the unit of [GeV-cm2-sec-sr]~^. The energy E is in the GeV unit, and the normalization 
factor is fixed to be A*",^ = 0.66 in the present analysis. By using these results, we 
estimate the quantities reported by PAMELA and Fermi-LAT. They are expressed as 



^prim _|_ ^sec.bkg 

e+ e+ 



prim 



+ ^ 



sec.bkg 



prim.bkg 



^prim 

e+ 



sec.bkg 



for PAMELA, 



for Fermi-LAT, 



respectively. In this estimation we should use the parameters which satisfy the constraints 
from the neutrino oscillation data, the LEV and the WMAP data. Such examples are 
shown in Eigs. 1 and 2. Although the ambiguity exists in the choice of each value of 



these parameters, it is absorbed into the assumed ip 



lifetime r^^^ . It can be justified by 



tuning the free parameter bi as found from eq. ([12]) o Thus, the lifetime r^^^ is treated 
free parameter in this analysis. 



^^It is interesting that the required value for bi can be consistent with the ones which explain the 
hierarchy of the coupling constants and the masses |13| . 
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Fig. 7 Flux of positron and electron predicted by the model for the data of PAMELA (the left frame) 
and Fermi-LAT (the right frame). Relevant parameters are fixed as Mi — 3.2 TeV and — 300 GeV. 
The flux is plotted for three typical cases such as (ce : : Cr, t^jvi) =(0:1:1, 4.0 x lO^^s), (0:1:0, 
6.7 X lO^^s) and (0:0:1, 3.3 x lO^^^s). 

The fluxes predicted by the model are plotted in each frame of Fig. 7 for some typical 
values of r^^^ and Cq, listed in Table 2. The data of PAMELA [17j and Fermi-LAT [18] 
are also plotted in the corresponding frame. We can fit the predicted flux in all the cases 
shown in Table 2 to the data of PAMELA well. However, the situation is different in the 
Fermi-LAT case. Although the predicted flux in the ) and (b) can be fitted to the 

observed data well, the case (c) can not be fitted to the data. The reason is that the 
positron produced from is softer than the one from /i"*". We find that the best fit is 
obtained in the case (b) where ip^-^ decays to fi^ dominantly. This case is also allowed 
from a view point of the constraint of diffuse gamma ray [12]. Since we aim to explain 
both anomalies in this analysis, we need to suppose a large mass for ipNi- It requires large 
neutrino Yukawa couplings as seen in the previous part. However, if we confine our study 
to explain the PAMELA anomaly only, rather light ipNi can also work well. In that case 
neutrino Yukawa couplings need not to be so large. 

The heavier DM component tpNi has also a radiative decay mode to the lightest neu- 
tralino x- Its one- loop diagram is shown in Fig. 5. This decay associates a characteristic 
gamma which can be found through the observation of the cosmic gamma rays. It has 
a line shape spectrum at the energy {Mf — m^)/2Mi which corresponds to the endpoint 
of the gamma ray spectrum generated through the processes such as the bremsstrahlung 
and the inverse Compton scattering associated to the tpNi decay and also through the 
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hadronization, fragmentation and decay of the final states. 
The width of this radiative decay is calculated as 

,3 



Ml 



i^r+i^r) 



where A and E are defined as 
A 



r2 2 2 



gMi (cos^ 0^1 {M^, ml, ml^,mj) + sin^ e^I{Ml 
+g*m^ (cos^ er,I{ml, Ml, m|, m^^) + sin^ ^^/(mj, M^, m|, m^J) 
i3 = Qm^ (cos^ 0,rjl{ml, Ml, mj, m^^) + sin^ OT^I{m\, Ml, mj, m^_)) , 
+e?*Mi (cos^ e^I{Ml, ml, ml^,m\) + sin^ ^^/(M^, mj, , m|)) 

g = {g'Z,,+gZ,^){M^,B*)[Y.hl,cle-'- 
In these formulas we use the definitions such a; 



Tf 2 2 2 2\ 



1 


2m^(m2 — 




1 




ml{ml - 




1 





r l a d 
-'1 2 ' 2 



2 2 



m. 



I ■ ■ "0 "(J "C 

-(2 — ^7, — — 7T 



J a, c 

-'I 2 ' 2 

2 \ / 2 2 2 

^ / mt mz, mi 



2 I 9 ' 2 ' 2 

.rn^rn^ m|\ _ /ml m^ m^ 

3l 2' 2' 2/ 3l 2' 2' 2 



(19) 



(20) 



Ji (ai,a2) 



1-02 



1 — ai 

h («!, ct2, as) = dx 
Jo 



log 



ai — a;2 



1 - a;2 
x(l — X 



h (ai,a2,a3) 



x(l - ai) + (a;2 - «3) 

a:;(l — x) 
x(l - ai) + (a2 - as) 



— ) log 

ai 



1 + 
1 + 



ai — a2 



a2 



a\ — a2 



ai(l — ai 



xai — a2 



x{l - ai) + {a2 - as' 

X — as 
x{\ - ai) + (a2 - as' 



log 
log 



a2 — 


a\x 




as - 


- X 




a2 — 


a\x 




as - 


- X 





(21) 



The contribution from the Higgsino component can be neglected since it is proportional 
to the lepton mass and then small enough. 

If we use these formulas, we can estimate the diffuse gamma flux generated by the 
'?/'7Vi decay. For example, we could predict the monochromatic gamma ray flux generated 



Although the function I {ml, 
can check that it is not singular. 



?,m^) may be considered singular at — for example, one 



17 




2400 2800 3200 3600 
Ml [GeV] 



Fig. 8 The decay width for — t' X7 as a function of the mass of ipNi- 



through this DM decay in the Milky Way halo as 

where an integral is done over the DM distribution along a line of sight. All astrophysical 
uncertainty is contained in this integral. The radiative decay width T^{il)j:q^ — )■ xi) is 
plotted as a function of the mass of ip^^ in Fig- 8. In this calculation we adopt the case 
(b) in Table 2 and use the following parameters: 

A = 7.74 X 10"^ mo = 400 GeV, M3 = 9000 GeV, - v^a = ^, 

tan/3 = 10, /i^ = 3600 GeV, S = 250 GeV, = 300 GeV. (23) 

These are the same ones used in the estimation of both the relic abundance of ip^i and the 
positron flux generated by the t/'tVi decay. This result suggests that it may be observed 
at the proposed Cherenkov Telescope Array in the future [43j. If this line shape gamma 
flux is observed, we can consider that it is a signature of the model for its peculiarity. 

3.3 Direct detection of the neutralino 

Direct detection of the DM can clarify the nature of DM [M]. Several experiments to 
search its elastic scattering with nuclei such as CDMSII, XENONIOO and XMASS are 
now under going or will start in near future. Some of these experiments have already 
constrain the models. For example, a vast region of the parameter space in the CMSSM 
has been excluded [S^. Thus, it is crucial to address the discriminative features of the 
present model, which could be expected to be found through these experiments. In 
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our model there are two components of the DM, ifj^^ and x- Since V'Afi does not have 
interactions with nuclei at tree level, the scattering cross section is heavily suppressed 
by the loop factor. Thus, it is difficult to detect it in these experiments. On the other 
hand, the lightest neutralino x can be scattered with nuclei at tree level since it has the 
same nature as the ordinary neutralinos in the MSSM. However, the constraint on the 
X mass and its scattering cross section imposed by the relic abundance can be different 
from the one in the MSSM as discussed in the previous part. Since the model has two 
DM components, the relic abundance constraint should be satisfied by both of these as 
shown in eq. (Q. Therefore, conditions for the parameters relevant to the direct search of 
X can be changed from the one in the MSSM, although the interactions of x with quarks 
are same as the MSSM neutralino. This could give a new possibility for the direct search 
experiments, which is not allowed in the MSSM case. 

The spin independent scattering cross section between the neutralino x aiid the nucleus 
with the atomic number Z and the mass number A is expressed as [3i 



4^ = ^[Z/p + (A-Z)/„]^ (24) 
vr 

where no momentum transfer is assumed. In case of large squark masses <^ rriq, 
the effective couplings of the neutralino x with the proton (/p) and the neutron (/„) are 
written as 

E ^-^ + ^fTG E -■ (25) 

q=u,d,s ^ q=c,b,t ^ 

where fq is the scalar four-point effective coupling constant whose concrete expression can 
be found in ^J. fjn^ represents the matrix element of nucleon defined by (A^|gg|A^) = 
fTqMn/mq and /tg is expressed as fxc = 1 — /t,- Although we use the values of 

q=u,d,s 

/t„ ^ given in [3l] , we adopt the smaller value of fr^ which is given in , 

/t„= 0.023, /t, = 0.034, /t. = 0.02, ioi N = n, 

/t„= 0.019, /t, = 0.041, /t. =0.02, ioi N = p. (26) 

In the numerical calculation, we treat and as free parameters in the allowed 
range shown in Figs. 3 and 4. Other relevant parameters are fixed to tan /3 = 10 and 
= + 50 GeV where is the slepton mass. Squark masses are assumed to be 
heavy enough. These parameters are those used in Figs. 3 and 4. 
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Fig. 9 The left frame shows the ahowed region in the {m^,af}) plane for the certain relic abundance of 
X- The right frame shows the detection rate of x in the XENON target for the x mass and the certain 
relic abundance of x- The MSSM corresponds to the upper edge of the region colored by blue. 



We use micrOMEGAs [35] in the analysis of the spin independent cross section with 
nucleon and the detection rate of the xo The left frame of Fig. 9 shows the region in 
the plane of and afj which is predicted by the model for various values of the relic 
abundance of x- The blue region corresponds to 0.10 < ^^^"^ — ^-^^ which includes the 
MSSM. On the other hand, the green region stands for the one with ^^^'^ — ^-^^ which 
can be consistent with the WMAP data in this model. The CDMSII and XENONIOO 
bounds are also plotted by a violet, red and black solid line in this frame respectively |45]. 
The region consistent with the relic abundance required for x, one of the DM components, 
is found to be much extended in comparison with the MSSM case. This occurs since the 
relic abundance constraint becomes much weaker than the MSSM case such as Q^h"^ < 
f^MSSM^^ — 0.11. The right frame shows the detection rate of x expected in the XENON 
target for each value of and the relic abundance of x- K suggests that the detection 
rate can be decreased by order one compared with the MSSM long as h"^ and 

VLy^h? are comparable. These features could allow us to distinguish this model from the 
MSSM. 
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The code includes the contributions from the x-gluon interaction via heavy quark loops. 
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4 Conclusion 



We have studied the nature of the DM sector in a supersymmetric extension of the ra- 
diative neutrino mass modeL An anomalous U(l) symmetry is introduced to explain 
the hierarchical structure of the couphng constants and mass scales in the model. The 
spontaneous breaking of this symmetry can induce a new Z2 symmetry which guaran- 
tees the stabihty of the hghtest odd parity particle. As a result, the model has two DM 
components as long as R parity is assumed to be conserved. However, since one of these 
discrete symmetries which guarantee the stability of DM is not exact due to the anomaly, 
one DM component is unstable to decay through a hugely suppressed term which is non- 
perturbatively induced via the anomaly effect. These DM components could be detected 
through the indirect search of the yields of the decaying DM and the direct search of the 
elastic scattering from nuclei by taking account that the DM relic abundance is composed 
of these. Positrons generated by the decaying DM can explain the cosmic ray anomaly 
reported recently. Parameter regions predicted by the direct detection can be different 
from the MSSM case since two DM components may contribute the rehc abundance in 
the same order. If the line shape gamma is observed in the cosmic ray, we might confirm 
the model by combining it with the direct search of the DM. Forth coming experiments 
for DM can give fruitful information to the model. 
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